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ABSTRACT 

It has been suggested that the final stages of accretion in present-day giant elliptical 
galaxies occur in an advection-dominated mode. The poor radiative efficiency of this 
accretion solution accounts for the fact that massive black holes, which are believed to 
reside in the centres of these galaxies, do not have the luminosities expected if accretion 
from the hot interstellar medium occurs at the Bondi rate. We discuss the advection- 
dominated solution for the nucleus of M60 (NGC 4649) and show that accretion at 
the Bondi rate is consistent with the core flux from radio to X-ray. This solution 
allows for a black hole mass of ~ 10 9 M as required by independent arguments. The 
successful application of this model to both M60, and previously to the well known 
nucleus of M87, suggests that accretion of hot gas in an elliptical galaxy creates the 
ideal circumstances for advection-dominated accretion flows to operate. 
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1 INTRODUCTION 



Most nearby, large elliptical galaxies should host a massive 
black hole left over from an earlier quasar phase (Fabian & 
Canizares 1988; Fabian & Rees 1995). Giant elliptical galax- 
ies also contain an extensive hot halo which provides a min- 
imum fuelling level for accretion. This leads to a problem. 
Fabian & Canizares (1988) have determined a limit to the 
accretion rate for a sample of bright elliptical galaxies pro- 
duced using the classical Bondi (1952) formula. They con- 
cluded that if the radiative efficiency of accretion is ~ 10 per 
cent, then all such nuclei would appear much more luminous 
than observed. The limits on the luminosity of the galaxies 
mean that black holes masses are of at most a few times 
10 7 Mq contradicting quasar counts and integrated lumi- 
nosities which imply black hole masses above 10 s - 10 9 M Q . 
Low-luminosity radio activity is very common in ellipticals 
(Sadler, Jenkins & Kotanij 1989; Wrobel 1991, Slee et al. 
1994) indicating that indeed there is some kind of central 
engine in the cores of these galaxies. 

The inconsistencies described above, can be reconciled 
by considering a situation in which the intrinsic radiative 
efficiency of the accretion is low. This has been noted by 
Fabian & Rees (1995) in the light of recent discussions 
of optically thin advection-dominated accretion solutions 
(Narayan & Yi 1995a,b; Abramowicz et al. 1995, Narayan, 
Mahadevan & Yi 1995 and references therein) . In this mode 
of accretion, the accretion rate, M , is low, the ions remain 
at the virial temperature and the radiative efficiency of the 



low density accreting material can be very small. The energy 
released by viscous friction is advected into the black hole. 

Accretion of hot gas in an elliptical galaxy may create 
the ideal circumstances for a hot ion torus to form (Begel- 
man 1986) and thus an advection-dominated accretion flow 
(ADAF) to operate Q In particular the successful applica- 
tion of low efficiency accretion via an ADAF in the giant 
elliptical galaxy M87 (which contains a 3 x 10 9 Mq black 
hole - Ford et al. 1995; Harms et al. 1995 ) by Reynolds et 
al. (1997; Paper I) suggests that such mechanisms can be 
of wide importance. There, we found that the luminosity of 
M87 is consistent with the observations only if Bondi accre- 
tion from the hot gas involves an ADAF. X-ray observations 
provide us with core gas densities and therefore the means 
for estimating the Bondi accretion rate. 

Here we present another example of a giant elliptical 
galaxy, M60 (NGC 4649), in which we believe accretion 
could be taking place via an ADAF. The giant elliptical 
galaxy M60 (NGC 4649) is an optically luminous galaxy in 
the Virgo cluster. It possesses very weak jets and small ra- 
dio lobes indicating an active nucleus. The radio spectrum 
of the nucleus is weakly inverted (Fabbiano et al. 1987) as 
expected from an ADAF. 

A deprojection analysis of data from the ROSAT }iigh 
Resolution Imager (HRI) shows that the interstellar medium 



* We note that the gas in the central region of elliptical galaxies 
is very likely to share angular momentum from stellar mass loss 
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(ISM) in M60 has a central density n w 0.1 cm -3 for a sound 
speed c s = 300 km s -1 (C. B. Peres, private communica- 
tion). The resulting Bondi accretion rate onto the central 
black hole is then 

M Bondi = 3 x 10 3 ( ^f) ^ Meyr- 1 . (1) 

If we assume a black hole mass of ~ 10 9 and a radia- 
tive efficiency of 77 = 0.1, the Bondi rate for a radiative 
efficiency of 77 = 0.1 predicts a luminosity of the order of 
6 x 10 43 ergs -1 . We show here that the luminosity of the 
core of M60 does not exceed 10 41 ergs _1 and is probably 
less than 10 40 ergs _1 over all energies. 

In this Letter we present a detailed examination of the 
possibility that the massive black hole in M60 accretes via 
an ADAF. In particular, we compute the spectrum expected 
and show that it is consistent with observations for a likely 
black hole mass (» 10 9 M Q ) and reasonable mass accretion 
rates given by eqn.(l). In contrast to M87 (Paper I) the ob- 
servations of M60 are not highly contaminated by jet emis- 
sion and therefore provide tighter constraints on the model 
spectrum. 

In Section 2 we compile data from the literature on the 
full-band spectrum of the core of M60 and present some 
additional data on the X-ray flux from the core. Section 3 
describes some details of our ADAF model spectrum calcu- 
lation. In section 4, we compare the model spectrum with 
the data and find that accretion rates comparable with the 
Bondi rate allow a black mass of ~ 10 9 M at least. 

2 THE SPECTRUM OF THE CORE EMISSION 
2.1 The M60 data 

In order to examine the nature of the accretion flow in M60, 
we have compiled the best observational limits on the full 
band spectrum of the core emission. Our aim is to obtain 
good observational limits on the core flux over a wide range 
of frequencies rather than to compile a comprehensive list of 
all previous observations. Some contribution from the weak 
jets and in particular from the underlying galaxy are un- 
avoidable and so the derived spectrum should be considered 
an upper limit to that of the accretion flow at the core of 
M60. The data are summarized in Table 1. 

Very Large Array (VLA) imaging by Stanger & War- 
wick (1986) and Wrobel & Heeshen (1991) shows an ex- 
tended radio source dominated by a compact (< 4 arc-sec) 
core. The extended structure defines radio lobes powered by 
weak jets with linear dimensions ~ 2 kpc. The slowly-rising 
spectral nature of the high-frequency emission is well defined 
by the data of Fabbiano et al. (1987) obtained from the 100 
m radio telescope of the MPIfR at Bonn. 

Einstein observations of M60 reveal that the overall X- 
ray emission is extended with most of it arising from the hot 
gas at kT = IkeV (Stanger and Warwick 1986). 

We have examined the ROSAT HRI dataset in order 
to constrain the nuclear X-ray flux of M60. Here we present 
the ROSAT HRI data resulting from a 19 923 s exposure 
performed on 1995 June 21. Fig. 1 shows the ROSAT HRI 
image of M60. Consistently with the Einstein HRI observa- 
tion, the ROSAT HRI image shows that M60 is extended 
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Figure 1. The core regions of M60 (NGC 4649) as imaged in 
a 20ks exposure with the ROSAT HRI. The distribution of the 
X-ray emission in the central region is quite complex and asym- 
metric. The diffuse emission is from the hot interstellar medium. 
The image has been smoothed adaptively prior to contouring; the 
minimum number of counts over which smoothing occurred is 25. 
Contour levels are equally spaced on a logarithmic scale, start- 
ing at 8 X 10~ 2 count s _1 arcmin -2 and increasing by a factor 1.6 
between adjacent contours. 
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Figure 2. The contribution of a point source to the emission of 
M60. The solid line shows the PSF of the ROSAT HRI compared 
to the actual count rate per unit area. 



X-ray source with most of the emission coming from hot 
gas in the interstellar medium. The low luminosity radio 
activity though, indicates that there is some form of ac- 
tive central engine in the core of M60. In order to find an 
upper limit to the nuclear emission of M60, we have con- 
sidered an extended component modelled by a King pro- 
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Table 1. Summary of data for the core of M60. 



Frequency 
v (Hz) (10- 

1.4 X 10 9 
4.75 X 10 9 
1.07 x 10 10 

3.3 x 10 10 
5.45 x 10 14 

2.4 x 10 17 



vF v 



0.37 
1.14 
2.78 
10.5 
< 180 
<150 



reference 



Hummel ct al. (1983) 
Fabbiano ct al. (1987), Wrobel (1991) 
Fabbiano et al. (1987) 
Fabbiano et al. (1987) 
Byun ct al. (1996) 
this work 



notes 



Westerbork 
Effelsberg, VLA 
Effelsberg 
Effelsberg 
HST 
ROS "AT HRI 



file plus a point source modelled by the PSF of the HRI 
(see Fig. 2). Our upper limit (at 90 per cent confidence) 
predicts a count rate from the point source component of 
4.7 x 10 _3 cts _1 (Fig. 2). Assuming the spectrum to be a 
power-law with a canonical photon index F — 2.0 modified 
by the effects of Galactic absorption (with column density 
A^h = 3.0 x 10 20 cm~ 2 ); this count rate implies a flux den- 
sity at IkeV of F(lkeV) = 1.5 x 10~ 13 ergoni^s" 1 keV" 1 . 
This result is fairly insensitive to the choice of power-law 
index. 



3 ADVECTION-DOMINATED TORI 

It has been found (Rees 1982; Narayan & Yi 1995ab; 
Abramovicz et al. 1995) that the advection dominated so- 
lution can occur only at accretion rates below M cr i t ~ 
o?Mvaa, where a is the Shakura-Sunyaev viscosity prescrip- 
tion and MEdd corresponds to the standard Eddington accre- 
tion rate. At such low accretion rate the cooling timescale for 
the ions exceeds the inflow timescale and most of the thermal 
energy is carried into the event horizon. The plasma is much 
hotter than in the classical, radiatively-efficient, thin disk 
solution and the flow consists of a thick torus. We assume 
that electrons and protons are coupled only by two-body 
Coulomb interactions and that electrons radiate by synchro- 
cyclotron, bremsstrahlung and inverse Compton processes. 
Since the Coulomb coupling is very weak, the ions remain at 
the virial temperature and very little energy is transferred 
to the electrons. The disk, in this solution, has low radiative 
efficiency. 

For convenience, we rescale the radial co-ordinate and 
define r by 



R_ 

Its' 



(2) 



where 7? is the radial coordinate and Rs is the Schwarzschild 
radius of the hole and the accretion rate, 



M 



Mr 



(3) 



The ADAF model, which is based on that of Narayan & 
Yi (1995) is the same as that used in Paper I for M87 
(Reynolds et al. 1997). By assuming that the system is un- 
dergoing advection-dominated accretion, we can predict the 
radio-to-X-ray spectrum of the accretion flow. The amount 
of emission from the different processes and the shape of 
the spectrum can be determined as a function of the model 
variables: the viscosity parameter, a, the ratio of magnetic 



to total pressure, /3, the mass of the central black hole, M, 
and the accretion rate, rn. For the moment, we take a — 0.3 
and (3 = 0.5 (i.e. magnetic pressure in equipartition with 
gas pressure) . The electron temperature at a given point in 
the ADAF, T e , can then be determined self-consistently for 
a given m and M by balancing the heating of the electrons 
by the ions against the various radiative cooling mechanisms 
(within r < 1000, it is found that T c is approximately con- 
stant in the range of parameter space of interest and it is 
~ 2 — 3 x 10 9 K). To determine the observed spectrum, we 
integrate the emission over the volume and take account of 
self-absorption effects. We have taken the inner radius of the 
disk to correspond with the innermost stable orbit around 
a Schwarzschild black hole, r in = 3, and the outer radius 
to be r out = 10 3 (for which the characteristic 2-temperature 
ADAF solution is mantained). 

In Fig. 3 we show the spectrum of the advection- 
dominated disk for m = 6.3 x 10 8 , 10 9 , 1.6 x 10 9 where m 
is defined as m = Mj M and dotm given by the Bondi 
accretion rate in eqn. (1). The peak in the radio band is 
due to synchro-cyclotron emission by the electrons in the 
magnetic field of the plasma. The X-ray peak is due to 
thermal bremsstrahlung. The power-law emission extend- 
ing through the optical band is due to Comptonization of 
the synchro-cyclotron emission: more detailed calculations 
show this emission to be composed of individual peaks corre- 
sponding to different orders of Compton scattering. The po- 
sitions at which the synchrotron and bremsstrahlung peaks 
occur and their relative heights depend on the parameters of 
the model. The synchrotron radiation is self-absorbed and 
locally gives a Rayleigh Jeans spectrum, up to a critical fre- 
quency, v c . The bremsstrahlung peak occurs at the thermal 
frequency v ~ ksT e /h. 



4 THE DATA VERSUS THE ADAF MODEL 

Fig. 3 shows the data on M60 with predicted spectra for 
different black hole masses and the relative Bondi accretion 
rates. A strong feature of an advection-dominated solution 
is that the whole spectrum can be determined by few param- 
eters. With a and ft given the only two parameters left are m 
and m. In Paper I we had an independent determination for 
the black hole mass, m, of M87. The only parameter to be 
determined was therefore the accretion rate. A comparison 
with the observational limits confirmed that the physically 
consistent accretion rate for such a system was indeed the 
Bondi rate. Here we therefore use the Bondi formula and 
vary the mass m. We find that the radio bounds on Fig. 3 
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Figure 3. Spectrum of M60 calculated with an advection-dominatcd flow extending from r m i n = 3 to r max = 1 000. The parameters are 
a = 0.3, ft = 0.5. Three model spectra are shown: for (i) m = 10 8 ' 8 — solid line, (ii) m = 10 9 ' 9 — dashed line, (iii) m = 10 9 ' 2 —dot-dashed 
line and their corresponding m given by Bondi accretion formula, eqn.(l). The squares represent the radio measurements and the arrows 
the upper bounds for the optical and X-ray core emission of M60 taken from the different references explained in the text. For M60, a 
distance of 15.8 Mpc is assumed. 



allow black hole masses m w 10 9 , similar to that required by 
independent arguments for dead quasars. For such masses, 
the low luminosity of the core of M60 is consistent with the 
data at different wavelengths only if the accretion involves 
an ADAF. Interestingly, the radio data points reproduce the 
slope (oc u 2 ^) of the self-absorbed synchrotron radiation 
(Mahadevan 1996) very well. Such agreement is even more 
evident than in M87 (Paper I) where contributions from the 
powerful jet had to be taken into account. It is important to 
point out that radio observations (Slee et al. 1994) of many 
bright elliptical galaxies typically show rising spectra with 
an average radio spectral index of 1/3. 



5 SUMMARY 

We have shown that if M60 consists of an advection- 
dominated system accreting from the hot interstellar 
medium at a rate determined by the Bondi formula, then 
the upper limit on nuclear black hole mass is much higher, 
M ~ 10 9 Mq, than that determined by Fabian & Canizares 
(1988), M < fewxlO 7 M . 

The poor efficiency of this advective accretion mode ac- 
counts for the low luminosity observed and allows for black 
hole masses of the order of 10 9 Mq consistent with the idea 
that bright ellipticals like M60, host dead quasars. 

This conclusion is strengthened by the previous appli- 
cation of a similar model to the well known active nucleus 
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Figure 4. Spectra of M60 calculated for m = 10 8 ' 8 . The effects 
of varying a and j3 are shown, (i) the dashed line is obtained for 
a = 0.3 and /3 = 0.95, (ii) the solid line for a = 0.5 and (3 = 0.5 
and (iii) the dot-dashed line for a = 0.1 and fi = 0.5 



in M87, which has a nuclear black hole mass of 3 x 10 9 M 
and luminosity several orders of magnitude fainter than that 
expected from the Bondi solution (Paper I) . 

An important test for this model would be high fre- 
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quency radio/submm observations. These should reveal an 
abrupt spectral turnover above the inner self-absorbed fre- 
quency (~ 3x 10 11 Hz) predicted by the model. Observations 
of the turnover would constrain the allowable viscosity and 
magnetic parameters a and f3. Fig. 4 illustrates how different 
combinations of a and /3 are either ruled out by the available 
data or constrained by the higher frequency observations. In 
particular, we suggest that a and f3 should be anticorrelated 
for a given mass in order to have agreement with observa- 
tions. This means that for relatively low values of a (~ 0.1 
within the context of ADAF) require magnetic fields below 
their equipartition value (i.e. f3 > 0.5), for a black hole mass 
- 10 9 M . 

High radio frequency polarization and interferometry 
studies could indicate and resolve the torus perpendicular 
to the jet. Better constraints on the X-ray flux will only be 
obtained with the improved resolution of AXAF. 
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